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I. INTRODUCTION 



The accurate evaluation of energy levels in atoms and ions requires inclusion of both the elec- 
tron correlations and quantum electrodynamic (QED) effects. At present, in spite of significant 
theoretical effort [1] there is no universal computational method which treats accurately corre- 
lations and QED effects all together. For heavy few electron systems, electron interactions can 
be treated perturbatively within systematic QED approach |12|]. For light systems involving few 
electrons, the most fundamental approach is based on the expansion of energy levels in the fine 
structure constant a, and also in the electron-nucleus mass ratio. At the same time, electron corre- 
lations are treated accurately by the use of explicitly correlated basis sets. This method has been 
advanced significantly in the last years by the calculation of m corrections to helium energy 
levels pI], ma'' to helium fine structure [I4L 51], and ma^ corrections in three- and four-electron 
systems la, [Zilsl]. The achieved numerical precision is sufficient to determine nuclear properties 
from isotope shift measurements of transition energies. For example, the most accurate determi- 
nation of the deuteron charge radius comes from the measurement of the 1S-2S transition in hy 



drogen and deuteriurn [bl]. Recently a series of measurements of isotope shifts in helium 111 



nil 



and lithium isotopes [|l2l . 11311 together with the intensive calculations of theoretical energy levels 
[Q, 1^ is], [3] brought the accurate values of charge radii of short-lived nuclei with respect to the 
stable isotope. However, at present, theoretical methods are not capable of predicting energy levels 
with such an accuracy, which would make available the absolute determination of nuclear charge 
radii, except for the hydrogenic s yste ms, where accuracy is limited only by small ma^ (Z a)^ 
higher order two-loop corrections iIITII . 

The most important in the accurate calculation of (light) atomic energy levels is the precise rep- 
resentation of the nonrelativistic wave function. The frequently used explicitly correlated Gaussian 
(ECG) functions give very accurate nonrelativistic energies, but the wave function does not sat- 
isfy the cusp condition, and for this reason this representation cannot be used for the calculation 
of m and higher order corrections. Also, the estimation of numerical uncertainties within the 
ECG method is quite problematic. More difficult to use is the Hylleraas basis set [|18 . HQ], but 
its achieved accuracy exceeds significantly results with ECG for three electron systems. In order 
to solve accurately Schrodinger equation we use the Hylleraas basis set with the number of func- 
tions of about 10 000 for S-states and 14 000 for P-states. The calculations of matrix elements of 
nonrelativistic Hamiltonian, as well as relativistic and QED operators are performed analytically 
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using newly developed recursion relations for the Hylleraas integrals l2o[ |21. 

In this work we present the most accurate calculations of nonrelativistic, leading relativistic 
and QED contributions including finite nuclear mass corrections and nuclear polarizability, to 
energy levels of 2^51/2, 3^5*1/2, and 2^Pi/2 states of Li and Be+. As a main result, we obtain 
ionization energies of Li and Be+, transition energies 2^5*1/2 — 3^5*1/2 and 2^S'i/2 — 2^Pi/2, and 
corresponding isotope shifts. In comparison to the former work of Yan et al. Ol] our results are 
in agreement for transition frequencies, but in slight disagreement for the isotope shifts in Be+. 
For the Li isotope shift both Ref. [7] and this work are in agreement with jl^. In all cases our 



numerical precision is about an order of magnitude higher, with the tota 



higher order terms. The comparison to experimental values 1124 , 



uncertainty dominated by 



2711 . apart from agreement 



with experimental transition energies in both Li and Be+, reveals small discrepancy for the Be+ 
ionization energy with the NIST data MM- Together with these benchmarking results, from our 
theoretical predictions and measured isotope shifts [Q,[l3, 28], we obtain improved nuclear charge 



radii for various isotopes including halo nuclei Li and Be 



II. RELATIVISTIC AND QED CORRECTIONS WITH THE NUCLEAR CHARGE RADIUS 

The energy level E{a, r]) as a function of a and rj = —fx/mj^j = —m/ (m + mjy) is expanded 
in power series in its arguments 

E{a, T]) = ma' [^(^.o) + ^ ^(2,1) + ^2 ^(2,2)] ^ ^ ^4 ^^(4,o) ^ ^ ^(4,i)] 

' +r/£(5'l)] +m«6 [£(6,0) ^^^(6,1)] +^«7£(7,0)^ 

and each coefficient is calculated separately from the expectation value of the corresponding 
Hamiltonian or the operator. The leading terms £^'''^\ and E^'^''^^ result from nonrelativistic 

Hamiltonian for three electrons and the nucleus with Coulomb interactions between them, 

^ m rriN ^ ^ Tab 

a a a>b 

In order to calculate relativistic and nuclear recoil corrections for atomic systems with the finite 
size nucleus, one should at first properly define the nuclear charge radius. Therefore, let us con- 
sider the interaction of the particle having spin s and the charge e with the electric field. The 
leading interaction related to the finite size is 



SH=--[ (r4) 5'^ + is' Qe ] d^E\ (3) 



where (r^j^) is the averaged square of the charge radius, is related to the electric quadrupole 
moment Qhy Qe = {sn (2 sn — 1)) Q, and 



-s + - s s- 

2 2 3 



(4) 



For the point spin s = 1/2 particle (r^j^) does not vanish and is equal to 3/(4 m^) , it is the so-called 
Darwin term and it depends on the value of the spin s uM- This term is the source of ambiguity in 
the definition of nuclear charge radii ISOfl. If we want the point particle to have the vanishing charge 
radius, the Darwin term should be excluded from (r^j^) . However, the value of the Darwin term 
for the arbitrary spin point particle with the arbitrary magnetic moment is unknown. Therefore, 
we propose the universal definition of (r^j^) by Eq. dS]), and thus include the Darwin term within 
the charge radius. A similar, but even more complicated problem appears when QED effects are 
being included. Now, assuming the definition of charge radius by Eq. ([3]), the atomic Hamiltonian 
including relativistic corrections and neglecting magnetic moment anomaly for electrons is of the 
form Jail] 



HbP = ^ Ha + ^ Hab + Hn + ^ HaN^ 
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In practice small relativistic terms involving nuclear mass are treated perturbatively and in this 
work we neglect all relativistic 0(m/M)^ corrections. These terms become much more important 
in muonic atoms and cannot be neglected there. 
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Let us now consider leading QED corrections of order m a^, which also include the inelastic 
contribution Epo\ due to the nuclear polarizability. Since we do not consider the hyperfine struc- 



ture, the spin of atomic nucleus can be neglected, and the QED correction takes the form [|32, 
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5^(rafe) -—ma^P 

TT 



(12) 



where 



In kn = 



27:aZ{J2J^{r, 



{cf>\P 



lim 



dV0*(r) 



— 9(r - a) + 4 vr 5^{r) (7 + In a) 



(13) 



^(r). (14) 



The electron-electron terms have been simplified in the above, since 5^{rab) does not vanish only 
for singlet states, therefore Sa ■ Sb6^{rab) = —3/A6^{rab)- We included in HaN the leading log- 
arithmic contribution that comes from the nuclear self-energy, but neglected all nonlogarithmic 
[m/MY terms which are proportional to 5^{ra)- They are not known for a general nucleus since 
they depend on the nuclear spin and charge distribution within the nucleus. For the calculation 
of the isotope shift in Li and Be+ we again neglect all 0{m/MY terms. These terms are im- 
portant for muonic atoms, but their calculation requires proper definition of the nuclear charge 
radius including QED effects, and its relation to the charge radius obtained from a different type 
of measurements such as the elastic electron scattering off nuclei. 

The last term in Eg. (fTOl) . E^ox is the nuclear polarizability correction. It is significant for halo 
nuclei such as ^^Li 1II6I1 . or whenever the isotope shift transition in the optical range reaches sub 
MHz precision, for example in the 1S-2S transition in deuterium The nuclear polarizabil- 
ity correction is expected to be significant also for ^^Be, as this nucleus has the largest known 
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B{E1) line strength among all nuclei. For this reason, we calculate it using experimental 11341] and 
theoretical [13 511 data for the electric dipole excitation of the ^^Be nucleus. 

Considering m corrections, they are well known for the hydrogen. Results for few-electron 



atoms are expressed in terms of the effective Hamiltonian ll36ll . Corresponding calculations have 
been performed only for low-lying states of helium [Q]. Calculations for three-electron systems 
are at present too difficult and therefore we use an approximate formula on the basis of hydrogenic 
values. 
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(15) 



It includes dominating electron-nucleus one-loop radiative, radiative recoil and pure recoil cor- 
rections [37]. We neglect electron-electron radiative corrections and the purely relativistic correc- 
tions, as we expect them to be relatively small, of order 10%. The relativistic m corrections are 
also very difficult to calculate. Its neglect is the leading source of uncertainty in the theoretical 
predictions for transition frequencies. Similarly, the nuclear recoil correction m? /M a^, which is 
still significant for the isotope shift, is also estimated on the basis of the known hydrogenic value 
in the above formula. This introduced some uncertainty in the determination of the isotope shift, 
from which nuclear charge radii are obtained. 

Due to numerical importance, one calculates approximately m a' contribution which is known 
exactly only for hydrogenic systems [|37ll . 
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It includes one-, two-, and three-loop corrections, and values of A, B, and C coefficients may be 
found in [37]. Following Ref. 1321] these hydrogenic values of order ma^ are extrapolated to 
lithium, according to 



£'^'\Z) = [2S^''\lS,Z)+S^''\nX,Z -2 



(17) 



2(53(r))is,z + (53(r))„x,z-2' 
for X = S, and for states with higher angular momenta S^'^\nX, Z) is neglected. We expect 
this approximate formula to be accurate to 25%. This completes QED corrections to transition 
frequencies and the isotope shifts in light atomic systems. 



III. COMPUTATIONAL METHOD AND NUMERICAL RESULTS 



In the construction of the wave function we closely follow the works of Yan and Drake in nl% . 
The global wave function ^ for both S and P states is expressed as a linear combination of tp, the 
antisymmetrized product A of the spatial function and the spin function x, 

ij = ^[0(ri,r2,f3)x], (18) 

= ^[0i(fi,f2,f3)x], (19) 

0a(n,r2,r3) = 0(fi,f2,r3), (21) 
X = «(1) /3(2) a(3) - (3{1) a{2) a{3) , (22) 

with rii being non-negative integers, Wi E R+, and the subscript a = 1,2, 3. The matrix element of 
the nonrelativistic Hamiltonian Hq in Eq. ^ or of any spin independent operator can be expressed 
after eliminating spin variables as 

(V^I-f^olV'') = (2 0(ri,r2,r3) + 2 0(r2,ri,r3) - 0(r3,ri,r2) 

-0(r2,r3,ri) - 0(ri,r3,r2) - 0(r3, rs, ri)|iJo I^Vi. ^"2, ^3)) , (23) 
{iji\Ho\^',') = (2 0i(ri,r2,r3) + 2 0i(r2,ri,r3)-0i(r3,ri,r2) 

-C(^2,r3,ri) - 0^(ri,r3,r2) - 0^(r3, r2, ri)|i7o r2, r3)) . (24) 

While for S-states Hamiltonian matrix elements can be written in one form, for P-states with 
the help of an additional ri,r2,r3 permutation they can take two different forms: (03|ifo|03*) 
or (02l-f^o|03*)- Next, all these spatial matrix elements are expressed as linear combination of 
Hylleraas integrals, namely the integrals with respect to r j of the form 

f(n n n n n n ^ - f ^^^^ f ^^^3 n^,ri-W2r2-W3r3 

J[ni,n2,n3,n4,n5,ne) — / - — / - — / - — e 

J An J Att J Att 

r^2~l ^,"3-1 ^,"4-1 ^ns-l ne-1 /^cn 
'23 '31 '12 '1 '2 '3 ; V^-*^ 

with non-negative integers nj. They are performed analytically forni, 722, ^3 = 0, 1 and by recur- 



sion relations for larger rii using formulas derived in [|20n . These recursions give the most accurate 
numerical values of Hylleraas integrals among all the methods developed so far. Nevertheless, 
multiple precision arithmetics has to be used in generating the Hamiltonian matrix, in order to 
avoid near linear dependence of Hylleraas basis functions. 
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The total wave function is generated from all (p in Eq. (1231) with ni satisfying a condition 

6 

Y,n,<n, (26) 

i=l 

for VL between 3 and 12. For each Vt we minimize energy with respect to the free parameters Wi 
in Eq. (|23l) . In order to increase the accuracy of the nonrelativistic wave function, following Yan 
and Drake [|l9ll . we divide the whole basis set into five sectors (six sectors for P states), each one 
with its own set of Wi%. To avoid numerical instabilities, within each sector we drop the terms 
with Hi > ns (or < n^) and for 724 = drop terms with rii > n2 (or rii < 722). This division 
allows for a significant improvements of nonrelativistic energies by optimization of all, five for S 
and six for P states, sets of w/s. Numerical results for 2^51/2, 3^51/2, and 2^Pi/2 of Li and Be+ 
for various sizes of basis sets are presented in Table HI The results for the ground state of Li are in 
agreement with our previous evaluation [|22ll . Results denoted by 00 are obtained by extrapolation 
to the infinitely large (complete) basis set, by fitting the function X{^1) = Xo + Xi/Q,^ with some 
integer n. The similar fit is used for all other matrix elements, presented in the following tables. 

Calculation of relativistic corrections, which are given by Eqs. ©-dH), involve spin independent 
and spin dependent terms. The matrix element of spin independent terms are calculated according 
to Eq. (|23] ). while spin orbit terms are obtained for P1/2 state by using 



{i^al'^Qc- ffc\i^b)j=i/2 = i(^(j)a{ri,r2,r3) -2Q3X [0fe(ri,r2,r3) +0fe(r2,ri,r3)] 

c=l 

+ {Qi -Q2 + Qs) X [0fe(r2, ra, ri) + ^^(rg, r2, ri)] 

+ + Q2 + 4) X [06(ri,r3,r2) + 0b(r3,ri,r2)]y (27) 

and the result for P3/2 is equal to —1/2 of that for Pi/2- The tensor spin-spin interaction vanishes 
for both the P1/2 and P3/2 states. All these matrix elements include Hylleraas integrals with 
Hi = —1, which are difficult to obtain accurately. We use the one-dimensional integral form for 
/(—I, 0, 0; 0, 0, 0) and /(O, 0, 0; —1, n^, ne) and other f's with Ui = —1 are obtained by recursion 



relations 



2211 . Since these recursions are not stable numerically, we used quadruple, sextuple. 



and octuple precision arithmetics written by Korobov II38I] to avoid loss of the numerical precision. 
It was especially important for excited states. Individual results for various operators are presented 
in Table In] and the total relativistic correction in Table Hill By the symbol [. . .]mp in these tables, 
we denote the mass polarization correction, namely 



Jmp 



^ ^ a>b 
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TABLE I: Nonrelativistic energy and relativistic and finite nuclear mass corrections in Li and Be"'", = oo 
is a result of extrapolation 



n £:(2.0) £(2,1) £(2,2) ^(4,0) £-(4,1) 

Li2lSi/2 

10 -7.478 060323 7863 -7.779903 106 67 -1.801631491 -12.04991825 10.0109103 

11 -7.478 060323 8615 -7.779903 105 96 -1.801631553 -12.04991416 10.010935 1 

12 -7.478 060 323 889 7 -7.779903 104 98 -1.801631587 -12.04991345 10.0109404 
oo -7.478 060323 906(8) -7.779903 104 4(7) -1.80163162(4) -12.049913 0(4) 10.010945(4) 

Li3i5i/2 

10 -7.3540984210040 -7.646138 262557 -1.677971728 8 -11.871 1920 10.014596, 

11 -7.354098421 3021 -7.646138 262 527 -1.6779717580 -11.871 1779 10.014649 

12 -7.354098 421 379 9 -7.646138 262 612 -1.6779717757 -11.871 1715 10.014687 
oo -7.354098 421426(19) -7.646138 262 65 (3) -1.677 971789(11) -11.871 168(8) 10.01472(4) 

Li2lPi/2 

10 -7.410156 532150 2 -7.656 895 306 5 -1.806107 -11.8013710 9.685 503 

11 -7.410156532 5860 -7.656895 2176 -1.806275 -11.8013659 9.685230 

12 -7.410156 532 628 6 -7.656895 191 1 -1.806407 -11.8013635 9.685 356 
oo -7.410156532 665(14) -7.656895 176(9) -1.80651(8) -11.801362(2) 9.68543(8) 

Be+ 2i5i/2 

10 -14.3247631766163 -14.777 682 315 84 -3.6340569619 -43.6880393 26.932300 

11 -14.324763176725 -14.77768231448 -3.634057 065 9 -43.6880344 26.932382 

12 -14.324763176763 5 -14.777 682 313 66 -3.634057 095 3 -43.6880263 26.932370 
oo -14.324763176784(11) -14.777 682 313 0(5) -3.634057 110(11) -43.688023(4) 26.93237(2) 

Be+ 3I51/2 

10 -13.922 789268 385 7 -14.351840 883 53 -3.294071 80 -42.3443398 27.3080808 

11 -13.922 789268 5307 -14.351 840 888 66 -3.29407332 -42.3443238 27.3082248 

12 -13.922 789268 5542 -14.351840 890 68 -3.29407394 -42.3443184 27.3082106 
00 -13.922 789 268 570(10) -14.351 840891 8(8) -3.2940744(4) -42.344317(3) 27.30821(2) 

Be+ 2iPi/2 

10 -14.179 333292 3197 -14.345 507 657 -3.728 487 -42.123 997 4 24.352 898 

11 -14.179 333 293 202 9 -14.345 507 560 -3.728 679 -42.123 985 4 24.352 585 

12 -14.179 333 293 342 7 -14.345 507 534 -3.728 739 -42.123 9782 24.352 564 
00 -14.179 333 293 42(3) -14.345 507 52(3) -3.728 78(4) -42.123 976(6) 24.352 55(4) 



QED corrections include two new terms, and the Bethe logarithm, see Eqs. (fT3l) and (fT4l) . 
Matrix elements of the term involve Hylleraas integrals with = —2. Their calculation 
using recursion relations was presented in Ref. Iisl], and numerical results are presented in Table 
nil We note a slow numerical convergence of terms and its mass polarization correction. 
Bethe logarithms are far the most difficult in accurate numerical evaluation and in this work we 
use the results obtained by Yan et al. in [|7D. We note a weak dependence of Bethe logarithms on a 
state, all close to the IS hydrogenic value. 
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TABLE II: Relativistic and QED operators in Li and Be+. Implicit sum over a and sum over a ^ h pairs 



are assumed. Bethe lo| 


^arithms are that from Yan et al. 


17] 








operator 


Li 2Si/2 Li 351/2 


Li 2Pi/2 


Bo+ 25i/2 


Be+ 35i/2 


Be+ 2Pi/2 



PaPb 


0.301 842780 3(6) 


0.292 039 841 2(2) 


0.246 738 644(9) 


0.4529191361(7) 0.429 051623 1(7) 


0.16617423(2) 


PaPb/iE - HY'PcPd 


-1.499 788 83(4) 


-1.385 93196(2) 


-1.55977(7) 


-3.181 137 97(3) 


-2.865 022 8(3) 


-3.562 60(5) 


Pt 


628.449022(4) 


622.859 40(2) 


620.044 977(8) 


2146.52076(5) 


2108.368 02(9) 


2088.652 95(4) 


[Pa] mp 


41.176 62(5) 


40.628 83(7) 


40.628 87(6) 


115.923 1(2) 


112.589 8(3) 


154.322 2(4) 


V ^ab r-^^ J ^ b 


0.871 195 809(9) 


0.859 817 43(18) 


0.792 85159(4) 


1.819 804 86(12) 


1.774 364 56(14) 


1.22681728(11) 




mp 


J. 644 ozo Z.\p } 


o.UUz y 1 jyii) 


17.854 151(2) 


16.414 367(3) 


1 1 Tin 1 /I ^ /" 1 1 \ 




13.842 610 787(6) 


13.736 502 84(7) 


13.676 197 06(7) 


35.1U5 055 /i(o) 


T A CT7 Oil £.1 \ 

Wll 6(1) 


34.245 82096(5) 


[5^(r-a)]mp 


0.484589 3(6) 


0.487 894 7(6) 


0.647 914(5) 


Pt nT7 t^r\n '3//i\ 
U.Vz / oU / 3(4; 


U.yjz Izy (2) 


1.831450(7) 




0.544 324632 0(7) 


0.536168 418 9(4) 


0.532274 098 9(3) 


1.580 538 588(3) 


1.537 328 373(1) 


1.518 990086(6) 


[&^{rab)]rap 


-0.082 009 7(4) 


-0.078 406(2) 


-0.066725(2) 


-0.4j1 I 5L\L) 


A /inn n/'n\ 
-0.4U9 93d 9(9) 


-0.351752(3) 




87.276740 9(3) 


86.289 222 5(7) 


85.743 2612(7) 


222.628 511 9(4) 


^1-7 '7nnc7nc/"7\ 

21 /. 799 5 79 5(7) 


214.562 695 2(8) 




0.273413(5) 


0.198 05(8) 


0.28957(5) 


-7.5146(1) 


-7.579 45(7) 


-6.794 2(2) 


[lALlmp 


-1.627 46(5) 


-1.645 4(2) 


-1.892(5) 


-7.5146(1) 


-7.579 45(7) 


-6.794 2(2) 




-308.314 23(6) 


-305.939 0(2) 


-304.517 7(2) 


-910.9194(1) 


-897.084 8(3) 


-887.841 7(3) 


In(fco) - 21nZ 


2.98106(1) 


2.982 36(6) 


2.982 57(6) 


2.979 26(2) 


2.98162(1) 


2.982 27(6) 


[ln(fco)]nip 


-0.113 05(5) 


-0.110 5(3) 


-0.1112(5) 


-0.125 58(4) 


-0.1171(1) 


-0.1217(6) 


ra/r^ X Pa ■ CTa 






-0.125 946 352(50) 






-0.9691317(11) 


[ra/r^ X Pa ■ CTa]mp 






0.376388(3) 






3.043 394(15) 


"fab/rlf^ y-Pa-Oa 






-0.22464070(5) 






-1.659492 5(3) 


[fab/rlt^ X Pa ■ CTa]mp 






0.570585(4) 






4.532 68(9) 


"fab/rlf^ X Pb ■ CTa 






0.038 473 60(6) 






0.3608516(2) 


Kft/'^ab ^ P6 ■ o'olmp 






-0.213 52(3) 






-1.54982(10) 


'aAa X piv ■ CTa 






0.022 524 93(9) 






0.339008 2(2) 



IV. TRANSITION AND IONIZATION ENERGIES 

Numerical values of all coefficients for relevant transition energies are presented for lithium 
and beryllium in the Table Hill does not include there the nuclear polarizability contribution. 
It is a small correction, which results from modification of the nucleus by atomic electrons. The 
effect of nuclear polarizability in ^^Li was estimated in |l6| . and ^^Be+ is obtained here. While it 
does not affect the absolute transition frequencies much, it is significant for the isotope shift and 
the corresponding determination of the charge radii difference, see the next section. 

Obtained results for the energy expansion coefficients are in general in good agreement with the 
former calculation of Yan and Drake in ^ [isl, with corrections and improvements in [7]. We 
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have not yet confirmed their results for the Bethe logarithms IM and use them in our coefficients 
g(5fi) £(^''^\ We note that the present numerical precision of expansion coefficients is high 
enough, that the leading uncertainty of transition frequencies comes from £^^''^\ more precisely 
from the rough estimate of the relativistic (nonradiative) m correction, which is about 10% of 
the leading radiative correction. 

Using this expansion coefficients in Table Hill and atomic masses from Table |IVl one obtains 
transition frequencies which are compared in Table |V] to the previous calculations of Yan et al. 
[Q, [m , 15 ] and to the experimental results. Small differences with results of Drake and Yan are due 
to the better numerical accuracy of our results and the inclusion of the finite nuclear size correction 
in transition and ionization energies for Be+. In comparison to experimental values we observe an 
agreement for both the lithium atom and beryllium ion, with one exception. Namely, theoretical 
ionization energy is larger by 0.06 cm^^ than the experimental value, and this discrepancy was 
already pointed out in [[zl]. The relatively lower accuracy of theoretical results for the beryllium 
ion comes from the neglect of the nonradiative m correction, which significantly grows with 
Z. The direct calculation of this correction for the three electron system is a challenge. A simpler 
approach would rely on matching the high Z results for lithiumlike systems [45] with the low Z 
results obtained by the expansion in a. We note that in spite of the relatively large uncertainties 
coming from £^^'^^ the obtained results for transition frequencies are the most accurate so far, and 
no other approach allows one for the systematic calculation of all corrections in the low Z atomic 
systems. 



V. ISOTOPE SHIFT DETERMINATION OF NUCLEAR CHARGE RADII 

The isotope shift in atomic transitions come mainly from different masses of nuclei. Much 
smaller effect of order ~ 10^^ is due to different nuclear charge radii. Nevertheless, precision of 
isotope shift measurements and theoretical predictions is enough to derive charge radii from the 
comparison of the experimental and the theoretical data. Moreover, this determination of charge 
radii is far more accurate than the one obtained from electron scattering off nuclei and can be 
applied even to unstable nuclei such as ^^Li and ^^Be. From expansion coefficients for Li and 
Be"*" in the Table Hill one obtains contributions to the isotope shift. The example for ^^Li-^Li and 
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TABLE III: Expansion of energy in powers of a and i] in Li. The last column presents values in units cm 
for the ground state ionization energy for ''Li and ^Be+ with atomic masses from Table |IVl 8^^''^'^ is the 
finite size correction with rch(^Li) = 2.39(3) fm [39] and rch(^Be) = 2.519(12) fm [40]. 



energy 35i/2 - 25i/2 


2-Pi/2 — 251/2 


i.jT . iiai/2 


i.jT . iiai/2 [L,m J 


£■(2,0) 


0.123 961902 48(2) 0.067903 79124(2 


T i 

1 98 1 46 91 1 238n 2"! 


43 488 220 301 (6'l 


£•(2,1) 


0.133 7648418(7) 


0.123 007928(10) 


21 1 01 3 905 1 ((\\ 


— 3 621 708 


£(2,2) 


0.123 659 827(11) 


-0.00488(7) 


235 270 010(19'! 


000 316 


£(4,0) 


0.178 745(2) 


0.248 5514(10) 


n 240 540 1 (A^ 


2 81 1 269(5'! 


<?(4,0) 
^fs 


-0.666 695(r2j^) 


-1.045 608 (r^j^) 


— U.OlU '•^'^Vch/ 


n nnn *^on 
— u.uuu oyu 


£(4,1) 


0.003 78(4) 


-0.325 44(11) 


0.013 50(3) 


-0.000 012 


£(5,0) 


-2.193(3) 


-3.478(3) 


-2.860 8(6) 


-0.243 99(5) 


£(5,1) 


-1.46(2) 


1.30(3) 


-1.889(3) 


0.000 013 


£(6.0) 


-9.2(9) 


-14.4(1.4) 


-12.0(1.2) 


-0.0075(7) 


£(b,l) 


-38.(10) 


-46.(12) 


-49.(12) 


0.000 002 


£(7,0) 


115.(29) 


215.(54) 


147.(37) 


0.000 67(17) 


£•(2,0) 


Be+ 

0.401973 908 21(2) 0.145 429 883 35(5) 0.669196 938 370(19) 146 871.751363(4) 


£(2,1) 


0.425 8414212(13) 0.432174 796(11) 


0.701595 7713(4) 


-9.374 767 


£(2,2) 


0.339 982 7(3) 


-0.094 72(3) 


0.721945175(15) 


0.000 587 


>^(4,0) 


1.343 705(3) 


1.564047(3) 


1.892 557(3) 


22.11892(4) 


c-(4,0) 
^fs 


-4.416 477(1) (ri) 


-7.198 308(l)(r2j^) 


-5.937 280(l)(r2j^) 


-0.002 955 


£(4,1) 


0.375 84(2) 


-2.579 82(5) 


0.62019(3) 


-0.000441 


£(5,0) 


-13.168(2) 


-21.864(11) 


-17.742(4) 


-1.5131(3) 


£(5,1) 


3.19(2) 


30.92(13) 


28.75(8) 


-0.000149 


£(6,0) 


-81.(8) 


-132.(13) 


-109.(11) 


-0.068(7) 


^(6,1) 


-361.(90) 


-442.(111) 


-477.(119) 


0.000018(5) 


£(7,0) 


1036.(259) 


1787.(447) 


1360.(340) 


0.006 2(15) 



Be+-^Be+ is presented in Table |VIl We observe that the leading nonrelativistic contribution 
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TABLE IV: Atomic masses of Lithium and Beryllium isotopes and the atomic binding energy. 



Li isotope 


mass [u] 


Ref. 


Be isotope 


mass [u] 


Ref. 


6Li 


6.015122794(16) 




■^Be 


7.016 929 83(11) 


[42] 


^Li 


7.0160034256(45) 


[43J 


9Be 


9.01218220(43) 


[42] 


8Li 


8.02248624(12) 


[41J 


lOBe 


10.013 533 82(43) 


[42] 


9Li 


9.02679020(21) 


[41J 


iiBe 


11.02166155(63) 


[44J 


i^Li 


11.04372361(69) 




i^Be 


14.042 890(140) 


[42] 


Eu 


-7.281 au 






-14.669 au 





,-1 



|7],b-| 


25 


l,c-| 


26 


l,d-| 


24 


l,e-| 


271 





LP. 2Si/2 


351/2 — 25i/2 


2-Pl/2 — 25i/2 


^Li (this work) 


43487.159 0(8) 


27 206.093 7(6) 


14903.648 4(10) 


■^Li (the) 


43 487.158 3(10)" 


27206.093 0(10)" 


14903.647 9(10)" 


^Li (exp) 


43487.15940(18)^ 


27 206.094 20(10)'^ 


14903.648 130(14)-^ 


^Be"'"(this work) 


146 882918(7) 


88 231.919(5) 


31 928.734(8) 


9Be+ (the) 


146 882.923(5)" 


88 231.920(6)" 


31 928.738(5)" 


^Be"*" (exp) 


146 88286^ 


88 231.915^ 


31 928.744^ 



gives at least 99.9% of the total isotope shift. The relativistic recoil corrections are small but 
still important, while the theoretical uncertainty is dominated by rough estimation of S^^'^\ The 
nuclear polarizability correction is significant for both ^^Li and ^^Be"*", and its value is presented 
in Table rvTl Result from ^^Li was obtained in [llql . while for ^^Be^ it is calculated here using both 
experimental js^l and theoretical data [35] for the so called B{E1) function which is the reduced 
line strength for the nuclear El excitation, 



2 4:71 dB (El) 



(29) 



in units fm^ MeV ^. The kets |0iv) and \E) denote the ground state of the nucleus and the 
excited state with excitation energy E, respectively. The relevant formula relating electric dipole 
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nuclear transition moment with the shift of atomic energy levels is [|l6ll . 



(30) 



where Opoi is a weighted electric polarizability of the nucleus and is given by the following double 
integral 



Q^pol 



16 a 



1 



(K + 



dE-\{^^\d\E)\' 



1 + 



dw 



E 



w E'^ + w'^ 



1 



+ 



(31) 



{k + 1){k* + 1) \k + 1 k* + 1^ 
where k = a/1 + 2 im/w. The first integral over the nuclear excitation spectrum may involve a 
sum over discreet levels, as it is the case of ^^Be. This nucleus has an excited state with E = 0.320 
MeV and B{E1) = 0.116 fm^ and a continuum spectrum starts at Et = 0.504 MeV. The result 



for Opoi using the experimental n34\ or the theoretical data n35\\ is 

apoi = 39.7(40) fm^ = 6.90(69) 10"^ (32) 
and this value is used to obtain the shift of energy levels z/poi in Eq. (|5]) for Be+ in Table IVTl 

TABLE VI: Contributions to the ^^Li - ^Li isotope shift of 35'i/2 — 2S'i/2 transition, and to the ^^Be+ - 
^Be+ shift of 2Pi/2 ~ 2'S'i/2 transition, with excluding the finite size correction. The second uncertainty of 
Afthe is due to the atomic mass. 



correction 


Li(3Si/2 - 25i/2) [MHz] 


Be+(2Pi/2 - 25i/2) [MHz] 


Ai/(2,i) 


25104.520 2(1) 


31568.5773(8) 


Al/(2,2) 


-2.9679 


0.765 7(2) 


Az.(4-i) 


0.0378(4) 


-10.035 0(2) 


Ai/(5,i) 


-0.1064(15) 


0.8777(36) 


Ai/(6,i) 


-0.020(5) 


-0.092(23) 


Ai^poi 


0.039(4) 


0.208(21) 




25 101.502 8(64)(27) 


31560.302(31)(12) 


Ai^thc [V] 


25 101.470(22) 


31560.01(6) 



From the difference between experimental and theoretical isotope shift one determines nuclear 
charge radii by using 

At/exp - At/the = Cab (r^h a - r^hs) > (33) 
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with constant C obtained from Eq. Q with including logarithmic relativistic corrections to the 
wave function at the origin 



(34) 



Using isotope shifts as measured for Li in jl'A Il3|l and Be"*" in t2M, we obtain nuclear charge 
radii for corresponding isotopes in Tables IVIIl and IVinl Results for Li are slightly more accurate 
than our previous determination in |ll6(] due to the more accurate z/thc and nuclear masses. Our 
results for Be isotopes, agree with the recent determination presented in [28]. The uncertainty of 
our Az^thc comes mainly from 25% of Av^^'^^ and 10% of Az/poi. Nonnegligible are numerical 
uncertainties of Bethe logarithms and their mass polarization corrections. The uncertainty of C 
coefficients comes from the estimation of relativistic correction to the wave function at origin, 
which is about 25% of the logarithmic part. Nevertheless, uncertainties in dr"^^ come mainly from 
the experimental value for the isotope shift, and the uncertainty of the final comes mostly from 
the charge radius of the reference nucleus. As we have already mentioned, the direct determination 
of the charge radius from the absolute transition frequency is at present not possible for lithiumlike 
systems due to insufficient precision of theoretical predictions. 

TABLE VII: Summary of isotope shift determination of Li charge radii from 35i/2 — 25i/2 transition 
with respect to ^Li, r(''Li) = 2.39(3) fm 13911 . the first uncertainty of u^hc comes from unknown higher 
order terms, the second uncertainty is due to the atomic mass. 



isotope 


l/exp[MHz] [13] 


i/thc[MHz] 


C [MHz fm-2] 


<)>e\[fm2] 


rch[fm] 


6Li 


-11 453.983(20) 


-11 452.820 5(23)(2) 


-1.5719(16) 


0.740(13) 


2.540(28) 


^Li 


8 635.782(44) 


8 634.981 2(17)(9) 


-1.572 0(16) 


-0.509(28) 


2.281(32) 


9Li 


15 333.272(39) 


15 331.799 5(31)(12) 


-1.5721(16) 


-0.937(25) 


2.185(33) 


i^Li 


25 101.226(125) 


25 101.502 8(64)(27) 


-1.576 8(17) 


0.176(79) 


2.426(34) 
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TABLE VIII: Summary of isotope shift determination of Be+ charge radii from 2Px/2 ~ '^^1/2 transition 
with respect to ^Be+, r(^Be) = 2.519(12) fm |40], the first uncertainty of i^the comes from unknown 
higher order terms, the second uncertainty is due to the atomic mass. 



isotope 


i/exp[MHz] [28] 


z^thcLMHz] 


C [MHz fm-2] 


<^r-2j^[fm2] 


?^ch[fm] 




-49 236.81(88) 


-49 225.736(35)(9) 


-17.021(31) 


0.651(47) 


2.645(14) 


iOBe+ 


17 323.8(13) 


17310.437(13)(11) 


-17.027(31) 


-0.785(76) 


2.358(21) 


iiBe+ 


31 564.96(93) 


31560.302(31)(12) 


-17.020(31) 


-0.274(55) 


2.464(16) 



VI. SUMMARY 



We have calculated nonrelativistic, relativistic, and QED contributions to low-lying energy 
levels of Li and Be+ ions, including finite nuclear mass corrections. The computational method 
uses the Hylleraas basis set with the analytic integration technique. The obtained results are the 
most precise to date, with the accuracy limited mainly by higher order relativistic m corrections. 



Using the experimental results for the isotope shift of 3S'i/2 — 25*1/2 transition in Li 1112 , 



Ij and 



2-Pi/2 — 25*1/2 transition in Be+ [28], we obtain improved charge radii for Li and Be isotopes. 
We note the significance of the nuclear polarizability effect in ^^Li and ^^Be+ and relativistic 
correction to the wave function at origin for the determination of charge radii. 

The presented computational method is limited by unlcnown higher order relativistic and QED 
corrections, which become more significant for heavier nuclei. One possible solution for the 
charge radii determination for heavier nuclei is the spectrosopy of the four-electron ion, for which 
we thiiik, accurate calculations can be performed with the help of a Gaussian basis set with linear 



terms 146, 



47|] . Apart from nuclear charge radii, precise atomic spectroscopy may bring informa- 
tion about the magnetic moment distribution within nuclei. Indeed a measurement of the hyperfine 
splitting in ^^Be"*" [|48|] may give the size of neutron halo, which cannot be probed by other means. 
However, the interpretation of the shift of the hyperfine splitting in terms of the Bohr-Weiskopf 
effect is not obvious, due to possible large nuclear polarizability effects [|49|]. 

The significant advantage of the presented computational approach with Hyllerras functions is 
the ability to calculate higher order relativistic [|50D and QED corrections, although such a calcu- 
lation is not simple. We aim to obtain O(a^) and 0{a^) corrections to the hyperfine splitting in 
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order to investigate nuclear structure correction with halo nuclei, and also to verify the accuracy 
of simplified approaches such as the relativistic configuration interaction or the multiconfiguration 
Dirac-Fock method. 
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